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We describe the design and performance of IceTop, the air shower array on top of the IceCube neutrino detector. 
After the 2008/09 antarctic summer season both detectors are deployed at almost 3/4 of their design size. With 
the current IceTop 59 stations we can start the study of showers of energy well above lO^'^ eV. The paper also 
describes the first results from IceTop and our plans to study the cosmic ray composition using several different 
types of analysis. 



1. Introduction 

IceTop is the air shower array part of the Ice- 
Cube neutrino telescope [1] . It supports IceCube 
in terms of relative and absolute pointing and 
guard against high energy air showers being mis- 
reconstructed as ultrahigh high energy neutrinos. 
With nearly 1 km^ area and atmospheric depth 
of 690 g/cm^ IceTop is also an excellent tool for 
cosmic ray physics research. 

The combination of IceTop with the Inlce muon 
detector operates as huge three dimensional muon 
telescope that has big advantages in studies of 
the cosmic ray spectrum and composition. This 
device can 

• measure the cosmic ray energy spectrum by 
studies of the shower hits in the surface de- 
tectors 

• measure the cosmic ray spectrum and com- 
position by studies of the angular distribu- 
tion of the muon bundles in the in-ice array. 

• measure the spectrum and composition us- 
ing coincident IccTop/Inlce events 

IceCube can not measure the number of TeV 
muons under ice. It, however, measures the en- 
ergy released by these muons in their propagation 
in 1 km of deep ice. The average surface energy of 
vertical muons that reach IceCube with sufficient 
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energy to create signals is slightly higher than 500 
GeV. 

We briefly discuss the design and current status 
of IceTop, its calibration and performance, early 
air shower results, and outline the possible cosmic 
ray measurements with this detector. 

2. IceTop design 

IceTop will consist of 80 stations that are lo- 
cated at an average distance of 125 meters from 
each other close to the top of each IceCube string. 
Each IceTop station consists of two frozen water 
tanks. IceTop tanks are plastic tanks with white 
diffusive reflecting inner surface of radius 0.93 m 
and height of 1 m. The freezing is a long compli- 
cated process because the water in the tank has to 
be freed of the dissolved air which after freezing 
creates bubbles that make the tank volume not 
uniform in optical clarity. In spite of the aver- 
age antarctic summer temperature of -25°C, the 
freezing usually takes about two months. 

Each tank is equipped with two digital opti- 
cal modules (DOMs) identical to those used by 
IceCube in-ice array that collect the Cherenkov 
light emitted by the shower particles that enter 
the tank. Note that the shower gamma rays con- 
vert in the tank and the experiment measures the 
energy flow of the extensive air shower (EAS). 

A DOMs contain a 10" Hamamatsu R7081-2 
photomultiplier tube, onboard high voltage con- 
verter, flashers and onboard digitization clcctron- 
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ics: two 420 ns, 300 MHz chips, each with four 
channels with different gains [2J. The two tank 
PMTs are run on different gains, one on high 
(5x10^) gain (HG) and the other on low (10^) 
gain (LG) for better dynamic range. 

Events with signal above threshold only in one 
tank of the station are considered hits by indi- 
vidual cosmic ray electrons, gamma rays, muons 
or hadrons. Coincident hits in both tanks of a 
station are considered air shower events. 

3. Status of IceTop in 2008 

The deployment of IceTop tanks goes simulta- 
neously with that of Inlce strings. In the begin- 
ning of 2008 both IceTop and IceCube had 40 
stations and strings. Nineteen new stations and 
strings were deployed during the most recent sea- 
son and this year the experiment will be almost 
3/4 complete. 
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Figure 1. Map of the IceTop array. One sees 
tanks in stations deployed before the 2007/08 sea- 
son (crosses), during that season (x's) and in the 
2008/09 season, where only the planned positions 
of the stations are shown with circles. 

Figure [1] shows the plan of the array at the end 
of January 2009. The inter-station distances vary 



between 110 and 145 meters. 

3.1. IceTop triggers and filtering 

IceTop trigger requires coincidences between 
six DOMs in a time window of 5/iS. This is called 
Simple Muhiphcity Trigger (SMT). When this 
occurs the whole detector (all IceTop and Inlce 
DOMs) is read. The frequency of IceTop SMTs 
in 2008 was 17 Hz. The bandwidth for data trans- 
fer from South Pole is not large enough to send 
all triggered events to the North. For this rea- 
son filtering is required for the event fraction that 
is transmitted. The filtering also checks on the 
number of stations that participate in the trig- 
ger and requires at least 3 stations for STA-3 and 
8 stations for STA-8 events. The frequency of 
STA-3 events in 2008 was 12 Hz. AU coincident 
triggers between IceTop and Inlce SMT (STA-3 
+ InlceSMT) and 5% of the coincidence events 
of Inlce SMT with surface activity (Inlce SMT -f 
IccTop) are also transmitted. Table [T] shows the 
different filters that were active in 2008. 

Table 1 

Frequency and percentage of different trigger 
events transmitted from South Pole 



Event type 


frequency 


% transmitted 


STA-3 


12.0 Hz 


20 


STA-8 


0.5 Hz 


100 


STA-3 + Inlce SMT 


2.2 Hz 


100 


Inlce SMT + IceTop 


30.7 Hz 


5 



4. IceTop performance 

A very important component of understanding 
the experimental performance is the calibration 
of the detector, starting with comparison of the 
signals in high and low gain DOMs and the indi- 
vidual calibration of each DOM. DOMs are with 
different sensitivity that depends on features such 
as the total tank reflectivity, ice quality and tem- 
perature (that changes following the outside air 
temperature). For this reason IceTop had spe- 
cial muon calibration runs that did not require 
shower triggers. These runs allow us to compare 
the muon signals in each tank and convert the 
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charge of the DOM signal from photoelectrons to 
VEMs (Vertical Equivalent Muons) which makes 
the whole array response uniform. Fig. [5] shows 
the muon peak in one of the DOMs of station 30. 

Some tanks are studied with small muon tele- 
scopes set in different locations on top of the 
tank. Coincidences between the tank and the 
muon telescope are used to study the tank re- 
sponse to muons hitting the tank at different po- 
sitions. Five such histograms for DOM 61 of sta- 
tion 30 are also plotted in Fig. [2l 
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Figure 2. Response of a DOM to muon hits 
(shaded background) and to coincidences with a 
muon telescope at different locations on top of the 
tank. 

One VEM value for this DOM is around 170 
photoelectrons depending on the exact fit of the 
distribution. The peak of the distribution is 
slightly higher than those of most muon tele- 
scope coincidences as the tank is hit by inclined 
muons that have longer pathlength in the tank 
and thus generate more light. The use of VEM 
in the shower analysis and the muon calibration 
makes the detector uniform. A new procedure is 
currently implemented that will allow us to have 
continuous muon calibration. 

Another important calibration is that of HG 
and LG DOMs. It is controlled by comparing 
the signal strength frequency reported by HG and 
LG DOMs which in the ideal case show smoothly 
joining distribution like the one shown in Fig. [3) 



Figure 3. Frequency of different VEM signal 
strength reported by HG and LG DOMs. Low 
gain DOMs are responsible for signals exceeding 
10 VEM. 

The comparison of the HG and LG DOMs in 
the same tank is used for calibration of individ- 
ual pairs before the overall distribution shown in 
Fig. [3] is made. 

5. First results from IceTop 

The first results on the cosmic ray spec- 
trum were obtained by the cosmic ray working 
group are best presented in the PhD thesis of 
S. Klepser [4|5] from the 2007 IceTop array that 
consisted of 26 stations. This array was fully ef- 
ficient for proton showers with energy above 1 
PeV. 

The work started with a study of the lateral dis- 
tribution of the air shower signal in VEM, which 
was fitted from experimental data as 

where r is the perpendicular distance from the 
shower axis, Sref is the expected signal at dis- 
tance Rref and /3 is a slope parameter related to 
the shower age. From simulations k was found to 
be constant at 0.3. 

As shown in Fig. [4] this form is quite different 
from the electron density as described in Greisen's 
formula. The reason is that IceTop tanks mea- 
sure also shower gamma rays and muons that that 
have different lateral spreads from the shower 
electrons. 
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Figure 4. Lateral spread of the shower VEM 
signal from a proton shower in IceTop tanks 
compared to Greisen's electron and muon lateral 
spread for showers of energy 10^ GeV. 



Sref equal to the average grid of 125 m was 
used in the analysis to relate to the proton pri- 
mary energy Eq- A response matrix was calcu- 
lated to relate the primary and reconstructed en- 
ergy using 5*125 in the case of protons and iron 
showers. At energy of 3 PeV the accuracy in the 
shower core determination was 9 meters and the 
relative error in logi^{E) was 0.05 for zenith an- 
gles less than 30°. 

Showers in three zenith angle ranges were an- 
alyzed. Assuming pure proton or iron compo- 
sition led to different primary spectra in the 
three angular ranges, while the assumption of 
mixed composition led to a single spectrum with 
a knee at 3.1±0.3(stat.)±0.3(syst.) PeV and 
spectral indices of 2.71±0.07(stat.) below and 
3.11±0.01(stat.) above the knee with systematic 
uncertainty of 0.08. The absolute normalization 
of the spectrum is lower than in most other mea- 
surements. 

This research was conducted when the IceTop 
detector Monte Carlo code was in a very early 
stage of development. The current version of the 
code is much improved and currently a new anal- 
ysis is being performed. 



6. Studies of composition sensitive param- 
eters 

There are several composition sensitive param- 
eters that can be measured by IceCube as a 
cosmic-ray detector: 

• The energy deposited by the shower muon 
bundle in 1 km of ice. Note that at primary 
energy 10^^ eV the average muon bundle in 
iron showers deposits about 2.4 times more 
energy than that of proton showers. 

• The waveforms in IceTop tanks can be used 
for a study of the GeV muon density, at 
least far away from the shower axis. 

• The lateral distribution of the shower signal 
is flatter for showers of heavy nuclei. 

• The risetime of the air shower front depends 
on the depth of shower maximum Xmax^ 
and is different for proton and iron show- 
ers. 

• The angular distribution of the showers in 
IceTop depends on the cosmic-ray chemi- 
cal composition as the showers initiated by 
heavy nuclei are absorbed much faster as a 
function of zenith angle. 

• The angular distribution of muon bundles 
in-icc in non-coincident events is also sensi- 
tive to the composition since proton showers 
generate higher energy muons that can pen- 
etrate to the detector at high zenith angles 
compared to those of heavier primaries. 

Our strategy will be to use as many complemen- 
tary observations as possible in order to decrease 
some of the uncertainties inherent in the inter- 
action models needed for simulations and inter- 
pretation of the data. All the parameters listed 
above will be studied, and the first four, which are 
applicable to IceTop and Inlce coincident events, 
will be determined for each coincident event. One 
should not forget that coincidental events are very 
close to vertical. 

Angular dependence of showers in IceTop 

As demonstrated by the preliminary analysis [5] 



5 



the angular dependence of rates on the surface is 
sensitive to composition. We are in the process 
of analyzing the full run of the IceTop 26-station 
configuration in 2007, which contains sufficient 
statistics to extend the analysis above 10^^ eV. 
This analysis uses cuts, such as containment of 
the air shower core in the perimeter of the array. 

We will continue to use this method for the full 
IceTop. It has two important features, the larger 
acceptance as compared to coincident events and 
its different dependence on models of hadronic in- 
teractions. With the full kilometer squared array 
we expect sufficient statistics to cover the 1-2 EeV 
energy bin. The method depends on shower ab- 
sorption and is closely related to a measurement 
of depth of shower maximum. It is therefore com- 
plementary to the analysis with coincident events 
that depends on the ratio of the signal of high 
energy muons to shower size. 

Study of muon bundles in coincident 
events: The acceptance of the telescope formed 
by IceTop and the deep IceCube neutrino tele- 
scope is approximately 1 /3 km^sr, sufficient to de- 
tect tens of events events per year above 10^* eV. 
If Ref. [3] is correct, we may find evidence sugges- 
tive of a transition from galactic to extra-galactic 
cosmic rays below this energy. Although full ac- 
ceptance will only be achieved at the end of the 
deployment period, already after the 2008/09 sea- 
son deployment IceCube is large enough to ex- 
plore fully the knee region up to 10^^ eV. Ex- 
pected event rates per year (assuming 70% effi- 
ciency) are indicated in Table [5) 

Table 2 

Air shower rates in IceCube in 2008 (IT4G) and 
upon completion (IT80) 



E (PeV) 


Rates yr-Mn(E) 




IT40 IT80 
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2 X 10« 10^ 
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Figure 5. The total number of photoelectrons 
(NPE) in IceCube is plotted as a function of the 
reconstructed air shower energy in IT40. 



As an example of what could be achieved by 
further analysis we show in Fig. [5] a scatter plot 
of the total number of photoelectrons (NPE) In- 
lce versus the shower energy reconstructed from 
the a month of IT40 data. The black dots show 
the average NPE as a function of energy. The cuts 
in the reconstruction require shower containment 
on the surface and Inlcc plus a good reconstruc- 
tion quality. The NPE is not yet corrected for the 
distance between the DOMs and the bundle tra- 
jectory. With the statistics shown we can study 
the composition in the energy range of 10^'^ to 
10^'^ GeV. The highest energy events collected 
already approach 10^ GeV. 

The data and analysis underlying this figure 
also illustrate the power of coincident events for 
energy calibration in IceCube. Events in which 
the shower core is inside IceTop can be recon- 
structed and energy assigned with a statistical 
uncertainty of order 5E/E = 0.1. To the extent 
that the muon content of air showers in the en- 
ergy range from 1 to 100 PeV are well understood, 
these coincident events can be used to calibrate 
the in-ice signal expected for deposition of an 
amount of energy expected from the muon bundle 
in the deep array. The study of muon bundles In- 
lce is complicated by the very large fluctuations 
in the high energy muons energy loss. 
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For downward events, where the signal is dom- 
inated by atmospheric muons, one can use the 
muon intensity (to the extent that it is known) to 
estimate the effective attenuation length (a com- 
bination of scattering and absorption 6J). This is 
done simply by calculating the expected intensity 
of muons at each depth from a standard param- 
eterization 7 (accounting for propagation in the 
ice |8l9j ) and fitting the observed DOM counting 
rate at each depth with a parameter that repre- 
sents the attenuation length of Cherenkov light 
at that depth. We have checked that the effective 
attenuation length derived in this way agrees well 
with direct measurements :6j. Such a simplified 
analysis is useful and fast for simulation of down- 
ward muon bundles where essentially all events 
are signal. 

The detailed study of muon bundles may also 
help us identify the production of very high pT 
muons in high energy showers. 

Identifying surface low energy muons from 
the detected waveforms: In addition to our 
principal discriminator of primary composition, 
which is the ratio of shower size at the surface 
to muon signal in the ice, we also have the pos- 
sibility of recovering some information about the 
relative contribution of GeV muons to the signal 
at the surface, which is also sensitive to compo- 
sition. Fig. [6] shows the waveforms of the signals 
in an air shower of reconstructed energy 5x10^ 
GeV seen by a detector at distance of 50 m from 
the shower axis and one at a distance of 250 m. 

The reconstructed energy of this shower is 
5x10^ GeV and zenith angle of 40°. The wave- 
form from the tank that is closer to the axis con- 
tains 11.34 VEM and width consistent with the 
shower front thickness. The waveform at 250 me- 
ters from the axis contains 1.13 VEM and is much 
more ragged. It seems possible that a muon hits 
first the tank and few low energy electrons and 
photons hit the tank later. We will investigate 
such waveforms coming from Monte Carlo events 
when the IT40 Monte Carlo produces sufficient 
statistics. 

One can in principle deconvolve the waveforms 
and estimate the number of muon peaks con- 
tributing to the signals in the tanks. Such a 
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Figure 6. Two waveforms taken at distances of 
50 m (left, with a LG DOM) and 250 m from the 
shower core with a HG DOM. The reconstructed 
shower energy is 5x10^^ eV. Arrows show the 
leading and trailing edges of the waveforms. 

treatment will be most effective in large air show- 
ers and at a significant distance from the shower 
core where the relative contribution of muons to 
the signal is higher. Thus such an analysis will 
be most appropriate when a large fraction of Ice- 
Top is deployed. Close to the shower core, as in 
Fig. [6l the estimate of the muon number may be 
more difficult. The fine timing resolution of the 
IceCube DOMs (3-4 ns per bin) may allow us to 
identify the GeV muons in air showers of primary 
energy above 10^^ eV. 
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